AFCRC-TN-55-797 


it,' 
to 


>^.     -^^    T^.     NEW    YORK    UNIVERSITY   M^^^k 


JL  iV  IX    ^     Institute  of  Mothematical  Sciences^ 
>-^       <<      Division  of  Electromoanetic  Research 


^^VJ  "^      ^•.        iri^  ...»  f^rrrUTE  OF  MAT' ?n:M  A  TCAL  SCIENCES 

>.      v.-^       ,s      Division  of  Electromagnetic  Research  f>i 

^^cccxx**-  3 

IS  1^1 1,  ■■!!■  r'lwwy  New  Yofk  3,  N.  Y,  ^ 

RESEARCH    REPORT    No.    EM-84 


Diffraction  by  a  Circular  Aperture  at  High  Frequencies 


HAROLD  LEVINE 


CONTRACT    No.    AF    19(122)-42 
SEPTEMBER,    1955 


NEW  YORK  UNIVERSITY 
Institute  of  Mathematical  Sciences 
Division  of  Electromagnetic  Research 


Research  Report  No.   EM-84 


DIFFRACTION  BY  A  CIRCULAR  APERTURE  AT  fflGH  FREQUENCIES 


Harold  Levine 


fT«'-vy\l<Ju       LCvAtA.^ 


Harold  Levine 


/\ 


Morris  Kline 
Project  Director 


The  research  reported  in  this  document  has  been  made  possible 
through  support  and  sponsorship  extended  by  the  Air  Force 
Cambridge  Research  Center,   under  Contract  No.    AF  19(122)-42. 
It  is  published  for  technical  information  only  and  does  not  neces- 
sarily represent  recommendations  or  conclusions  of  the  sponsor- 
ing agency. 

New  York,    1955 

NEW  YORK  UIN^RSJTX 

INSTITUTE  OE  MATHEN4ATICAL  SCIENCES 
LIBRARY 

75  Waverly  Phot,  New  York  3.  N.  Y. 


-  i  - 


Abstract 


An  investigation  is  made,  based  on  integral  equations  and  variational  prin- 
ciples, for  the  purpose  of  determining  the  first  correction  to  the  plane  wave 
geometrical  transmission  cross-section  of  a  circular  aperture  in  an  infinite 
plane  screen.  The  correction  is  sensitive  to  the  screen  boundary  condition, 
and  in  this  paper  details  are  given, for  nonnal  incidence  only,  for  the  acoustically 
'soft*  screen,  where  the  wave  function  vanishes.  There  are  two  ways  of  proceeding; 
1)  to  use  equations  which  characterize  the  normal  derivative  of  the  wave  function 
at  the  screen,  and  2)  to  start  with  the  wave  fvmction  itself  in  the  apertiire.  As  a 
consequence  of  approximations  employed,  there  is  some  disagreement  in  the  results. 
The  screen  formulation  turns  out  to  possess  advantages,  and  predicts  a  cross-section 

/     sin(2ka-  J)\ 
\^  V5r    (ka)^/2  J 

where  a  is  the  aperture  radius.  In  either  approach,  use  is  made  of  the  well-known 
problem  of  diffraction  by  a  half -plane  or  straight  edge,  a  configuration  irtiich  is 
relevant  when  the  aperature  dimension  greatly  exceeds  the  wavelength.  It  is  shown 
that  the  straight-edge  screen  distribution  due  to  a  grating  plane  wave,  and  not  the 
primary  wave,  plays  the  decisive  role  for  the  correction  in  the  cross-section.  A 
comparison  with  other  approximation  techniques  is  given,  and  emphasis  is  placed  on 
the  role  of  boundary  conditions  at  the  aperture  rim. 
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!•   Introduction 

The  mathematical  analysis  of  diffraction  phenomena  has  been  significantly 
extended  in  recent  year«.     A  variety  of  methods,  new  and  old,  have  proved  their 
worth  in  enlarging  the  roster  of  problems  amenable  to  calculation,  particularly 
as  regards  configurations  of  geometrical  complexity.     It  should  be  remarked, 
however,  that  the  form  of  solutions  thus  obtained  may  not  always  be  well  suited 
for  describing  certain  features  of  the  scattering  problem,   such  as  the  high- 
frequency  or  short-wavelength  characteristics.     The  latter  can  be  recognized  in 
the  zero-wavelength  (or  geometrical  optics)  limit,  although  without  the  inter- 
esting distinction  that  is  manifest  at  long  wavelengths,   say,  between  like 
scatterers  of  different  composition.     To  exhibit  a  corresponding  difference  in 
the  scattering  at  high  frequencies  one  must  refine  the  geometrical  optics  approxi- 
mation. 

By  vay  of  illustration,  let  us  consider  plane  soxond  waves  (of  harmonic  time 
dependence  and  infinitesimal  amplitude)  impinging  normally  on  a  circular  aperture 
.2 


Transmission  coeffcients  for  a  circular  aperture  for  nonnal  incidence. 
t]_     -     acoustically  'soft'    screen 

t2     -     acoustically   'hard*    screen 
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in  a  plane  screen.     In  Fig.  1,   the  energy  transmission  cross-section  of  the 
aperture   (relative  to  its  area)  is  depicted  for  both  an  acoustically  soft  (I) 
and  hard  (II)  screen.     If  the  acoustical  designations  are  interchanged,  these 
curves  also  yield  values  of  half  the  total  scattering  cross-section  of  a  cir- 
cular disk.     When  the  wavelength  A.  is  large  compared  to  the  apert\ur«  radius  a, 
there  is  an  evident  contrast  in  the  transmission  ci-css-sections  for  the  two 
screen  boundary  conditions.     The  origin  of  this  contrast  is  readily  traced  to 
the  different  orders  of  magnitude  for  the  incident  field  and  its  normal  deriva- 
tive, at  the  aperture  plane;  the  latter  is  smaller  by  the  amount  1/X.     It  is 
also  clear  that  the  boundary  condition  affects  the  manner  in  which  the  cross- 
sections  approach  the  comnon  geometrical  limit  (X  »  0),     The  exact  series  re- 
presentation of  each  cix)S5-section  >•  *   J ,  on  the  basis  of  which  of  the  curves 
in  Fig.  1  are  sketched,  lends  itself  to  approximation  at  long  wavelengths* 
However,  such  representations     prove  uns\iitable  at  short  wavelengths,  where 
even  verification  of  the  georietrical  optics  limit  is  analytically  diffic\ilt, 
^ry  further  comparison  of  the  two  boundary -value  problems  when  X.  »  a  is  thus 
far  unavailable. 

To  resolve  features  such  as  pointed  up  in  the  preceding  discussion  is  one 
aspect  of  a  study  concerned  with  the  behavior  of  solutions  of  the  wave  equation 

(1)  (V^  ♦  k^)  T(^   -     0 

for  ^  '  -r-  ->cxD.  The  nature  of  asymptotic  expansions  in  inverse  powers  of  k 
for  wave  functions  has  been  considered  by  Luneburg  i-  ■' ,  Kline'-  •',  and  more  re- 
cently by  Keller  and  Friedlander "^  K     These  authors  obtain  and  discuss  condi- 
tions for  successive  terms  in  the  asymptotic  expansion,  of  which  the  first  cor- 
responds to  geometrical  optics.  The  technique  is  applicable  in  some  scattering 
problems'-  -*,  although  th«»t  of  the  aperture  is  not  at  present  within  its  scope. 
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It  is  the  purpose  of  this  and  a  ccmpanion  paper  to  study  the  high-fr«quency 
transndssion  pixjperties  of  circular  apertures  in  soft  and  hard  screens,  with  the 
aim  of  clarifying  the  distinction  manifest  in  Fig.  1,  and  deducing  the  lowest- 
order  deviations  (in  reciprocal  powers  of  ka)  from  the  geometrical  cross-section. 
A  suitable  framework  for  such  considerations  is  provided  by  integral  equations  and 
variational  principles,  whose  fonns  depend  on  the  particular  boundary- value  prob- 
lem j  it  is  thus  appropriate  to  deal  separately  with  the  two  screen  types.     In  this 
paper  we  shall  be  concerned  with  the  analysis  of  transmission  through  an  aperture 
in  a  soft  screen.     The  problem  of  the  hard  screen  is  not  treated  here,  althovigh 
some  comments  on  it  are  included. 

Section  2  contains  preliminary  information  relevant  to  the  boundary-value 
problem;  various  high-frequency  approximation  procedures  are  described  and  the 
results  previously  obtained  are  given.     Section  3  is  devoted  to  an  integral  equa- 
tion characterizing  a  screen  distribution,   and  Section  U  to  its  conpanion  varia- 
tional principle.     The  aperture  integral  equation  is  then  considered  (Section  5), 
and  finsLLly  an  alternative  technique  for  the  screen  basis  is  presented  (Appendix 
1). 

It  may  be  expected  that  similar  procedures     are  appropriate  for  other  config- 
urations,  although  the  special  simplicity  of  plane  scatterers  will  become  apparent. 
It  should  be  noted  that  the  present  method  has  the  advantage  of  dispensing  with 
specific  functions  related  to  the  geometry  under  consideration  (spheroidal  functions 
for  the  circular  aperture'-  *   ■'), 


The  analogous  problem  involving  plane  electromagnetic  waves  and  a  perfectly 
conducting  screen  has  been  treated  by  H.  Chang  in  a  thesis  at  Harvard  University 
(1955). 
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2,       Description  of  methods  and  results 

Consider  the  infinitesimally  thin  screen  S,  located  in  Uie  plane  z  ■  0,  with 

ikz 
an  aperture  A  of  finite  extension.     The  scalar  plane  wave  e         is  assmned  to  be  in- 
cident on  the  aperture.     If  T(r)   denotes  the  acoustic  velocity  potential,  and  obey« 
the  wave  equation  (1),  we  have  T  ■  0  or  -r—  ■  0  on  the  screen  according  as  the  latter 
is  soft  or  hard.     In  terms  of  the  Gi^en's  function 


(2) 


G(^,i^') 


exp 


{ikft  -l^'l} 


Un    |r  -■?'! 


there  are  two  familiar  representations  of  T(r)  in  the  half -space  t  >  0: 


(3) 


T(r)  -  2 


1{P)   ^G(r,r') 


S+A 


z'  -  0 


dS', 


and 


(U) 


T(r)  -  -2 


az* 


T(?') 


'^S+A 


G(r,p')dS', 


where  r,  p  denote  points  that  lie  outside  and  on  the  apertxire  plane  respectively. 
Thus,  the  transmitted  wave  function  T(r)  is  determined  when  its  values  or  those 
of  its  normal  dei*ivative  are  specified  in  the  apertui*©  and  on  the  shadow  side  of 
the  screen.  Followine'  Braunbek^  -•,  we  may  tabulate  a  priori  our  knowledge  of  these 
boundary  values  for  a  soft  screen,  with  the  assumed  excitation: 

Table  I 


Boundary  values  of  f 
S 


Boundary  values  of  d  T 


T^act  0 

Kirchhoff  o 

approximation 


unknown 

1 


xinknown 
0 


Ik 
ik 
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Figure  2  (after  Braunbek) 


That  B  T  is  rigorously  equal  to  9_i^  in  A  is  a  consequence  of  symmetry 
in  the  geometrical  regions  (i.e.,  half-space)  on  opposite  sides  of  the  aperture, 
coupled  with  the  infinitesimal  thickness  of  S.  It  is  evident  that  the  Kirchhoff 
approximation  assumes  incident  or  unperturbed  distributions  in  the  aperture,  and 
if  used  with  the  representation  (U)  gives  rise  to  an  unspecified  error  in  tnin- 
cating  the  integral  beyond  A.  The  Kirchhoff  proposal  for  calculating  the  trans- 
mitted field  has  long  been  recognized  as  deficient  in  principle,  except  in  the 
geometrical  optics  limit.  When  the  wavelength  is  smsill  compared  with  aperture 

dimensions,  the  normal  derivative  8  Y  on  the  shadow  face  of  the  screen  is  signi- 

z 

ficant  principally  in  a  strip  region  bordering  the  aperture,  which  is  a  few 

wavelengths  in  width  (see  Fig,  2),  Braunbek 's  proposal  for  improving  on  the 

Kirchhoff  s^proximation  is  to  characterize  and  incorporate  in  (U)  values  of 

3  T  in  this  strip  which  pertain  to  a  local  Infinite  straight  edge  tangent  at 
z 

each  point  to  the  aperture  rim.     The  straight-edge,  or  half -plane  distribution 
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for  piane-wave  excitation  is  well  known  from  the  original  investigation  of 
Sommerfeld  l-  J,   and  the  work  of  many  subsequent  authors. 

In  a  similar  manner,  the  representation  (3),  which  involves  only  the  aper- 
ture values  of  T  when  the   screen  is  soft,   can  be  extended  beyond  its  Kirchhoff 
approximation.     To  do  this  one  must  take  into  account  an  additional  contribution 
to  the  integral  in  (3)  which  arises  from  a  strip  region  within  the  aperture  and 
near  the  rim,  where  the  wave  function  differs  significantly  from  its  incident 
value  J  indeed,   at  the  irLm  the   correct  wave  function  vanishes.     The  latter  differ- 
ence, not  known  exactly,  is  to  be  identified  locally  with  the  help  of  the  aper- 
ture distribution  appropriate  to  a  half -plane  whose  straight  edge  is  determined 
by  the  rim  tangent. 

For  the  aperture  of  circular  shape,  Braunbek'-  J  thus  obtains  the  following 
expressions  for  the  amplitude  of  the  transmitted  wave  en  the  axis  of  symmetry 
(the  z-axis),  for  the  representations   (U)  and  (3)   respectively: 

/®      ikr  /°°  ikr 

(5)  T  =  -ik  j       t^  p  dp  -  ik    I     F(^)  1;—  p  dp. 


r 
o  'a 


and 


(^)  '^  -  - 1  n  (V ) 


where 

(7)  i:  -  -  >|  -     k(p  -  a). 

and 


*  ikr  /"  ^        ikr 


pdp  ♦  j  $(7)  ~  (V-)P  dp. 


00 

2 


00  _ 

(9)  $(7)   -^     e-i^/^    j     e^^      dt. 
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^   Z 


Figure  ,3 


Here  the  aperture  radius  is  s,   and  the  geometrical  significance  of  r  and  p 
can  be  inferred  from  Fig.  3,     The  first  integrals  in  (5)   and  (6)  are  evidently 
Kirchhoff  terms,   although  the  former  actually  gives  the  exact  aperture  contribution 
(see  Table  I),     For  large  values  of  ka,  and  when  P  is  not  near  the  aperture  plane, 
the  functions  T  from  (5)  and  (6)  agree  as  regards  terras  of  lowest  order  in  recip- 
rocal powers  of  ka;  the  higher  terms  are  different  and  in  fact  outside  the  scope 
of  this  approximation.     In  loi*est  order,  then     (cf,    [9])^ 


(10) 


ikr 


T^e"^-e       °    /Up  , 


ka  »  1, 


where 


sin  ©. 


It  is  easy  to  verify  that  if  we  let  the  point  P  recede  from  the  apertiire   (z  ->  oo), 
then  the  expression  for  Y  goes  over  into  a  spherical  wave  form,  namely 
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ikz 
(il)  Y  i  A  — --     ,  z  ->oo 


where 

(12)  A ^ —    -  ^     , 

The  imaginary  part  of  the  amplitude  (12)  can  be  traced  to  the  Kirchhoff  integrals 
in  (5)   and  (6),  while   the  real  part  has  its  origin  in  the  straight-edge  corrections. 
Now,  if  we  use  the  well-known  theorem"-   -•: 

(13)  ^  ■  -  Y^     ^  ^» 

which  connects  the  plane-wave  transmission  cross-section  of  the  aperture  (for  normal 
incidence  here)  with  the  imaginary  part  of  the  forward  scattered  spherical  wave  ampli- 
tude, it  is  evident  that  in  the  approximation  (12)  we  have 

(111)  a-  =  na^  , 

the  geometrical  result.     Thus,  while  the  Braunbek  procedure  improves  the  scatter- 
ing amplitude  relative  to  the  Kirchhoff  theory,   the  relation  (13)  shows  that  the 
correction  is  not  sufficient  to  alter  the  cross-section.       The   cross-section 
can  naturally  be  identified  with  the  integrated  energy  flux  over  all  directions  in 
the  half-space  z  >  0;  however,  as  Braunbek  indicates,   there  is  some  difficulty  in 
obtaining  from  (5)  or  (6)   a  form  of  the  scattered  intensity  in  the  radiation  field 


*  From  the  highe/  terms  in  1/ka  that  the  representations   (5)  and  (6)  yield,  one 
can  obtain  frequency-dependent  additions  to   (lU)  via  (13),  but  these  lack  quan- 
titative  significance}  thus,  the  succeeding  term  is     0((ka)     ),  whereas  we  shall 


-3/2 
prove  that  it  should  actually  be     0(ka)  ), 
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which  holds  uniformly  for  all  directions,   and  consequently  he  gives  no  explicit 
correction  to   (lU). 

Since  the  cross-section  theorem  has  obvious  merit  for  computational  purposes, 
we  use  it  with  more  accurate  evaluations  of  the  forward  amplitude  A,  which  may  be 
accomplished  by  improving  on  the  Braunbek  distributions  in  the  plane  of  the   screen. 
The  lowest-order  correction  to  the  cross-section  (lU)  is  in  fact  obtainable  after 
extending  the  Braunbek  screen  distribution.       This  extension  is  also  derived  from 
a  straight-edge  configuration,  where  now  the  excitation  arises  from  a  secondary 
plane  wave  at  grazing  incidence.         The  grazing  wave  describes  an  interaction  be- 
t\«en  points  at  roughly  opposite  ends  of  an  aperture  diameterj  that  this  is  con- 
sistent with  a  correction  to  the  geometrical  cross-section  can  be  inferred  from  a 
ray  description,   according  to  which  undeviated  rays  passing  through  the  aperture 
provide  the  geometrical  part,  whereas  diffraction  effects  originate  with  rays 
scattered  from  the  rim  across  the  aperture  and  thence  along  the  incident  direction. 
Higher-order  corrections  evidently  refer  to  multiple  scattering  of  rays  across  the 
aperture . 

The  corresponding  description  of  the  aperture  distribution  is  rather  more 
difficult  to  obtain,   although  it  is  clear  that  the  straight-edge  variation  near 
the  rim,  as  generated  by  the  incident  wave,   does  not  provide  the  interference  or 
standing-wave  character  inherent  to  the  proper  aperture  distribution.     An  accurate 
form  of  the  aperture  distribution  at  short  wavelengths,  both  near  the   rim  and  at  the 
center,  would  seem  to  require  mo]*e  powerful  integral  equation  techniques. 

More  precisely,   there  is  also  a  geometrical  factor  that  modulates  the  straight- 
edge function,  with  maximum  value  unity  at  the  edge. 

S.  Karp  and  A.  Russek  have  considered  this  kind  of  interaction  in  the  case  of  an 
infinite  slit  formed  by  two  half -planes j  see  Research  Report  EM-75,  New  York 
UniTersity,   Division  of  Electromagnetic  Research, 
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Another  apn roach  to  the  detemination  of  the  scattered  amplitude  is  throxigh 
variational  princip'les  *-    J,     Here,  hovever,  we  are  confronted  by  the  formidable 
nature  of  the  variational  integrals  involved  when  the  Braunbek  or  more  refined 
distributions  are  utilized.     However,  the  integrals  can  be  evaluated  approximately 
when  the  Green's  function  occurring  therein  is  modified  in  a  manner  suggested  by 
the  derivation  of  such  distributions.     For  instance,  the  lowest-order  correction 
to  (ik)  emerges  from  a  variational  principle  with  use  of  only  the  Briaunbek  screen 
distribution,  if  the  Green's  function  is  suitably  characterized.     It  is  noteworthy 
that  even  a  first  correction  to  the  geometrical  cross-section  obtained  by  any 
method  proves  reliable   at  moderate  frequencies,  a  feature  welD  known  in  other  in- 
stances of  asymptotic  developments.     This  circumstance  contrasts  markedly  with 
the  behavior  at  low  frequencies,  where  successive   terms  in  a  power  series  expan- 
sion do  not  significantly  extend  the  qualitative  or  quantitative   aspects  of  the 
cross-section. 

The  explicit  correction  to  the  geometrical  optics  cross-section  (lii)   is  of 
interest  in  itself,   and  also  as  a  basis  for  comparison  between  different  methods 
of  approach.     Our  calculations,  given  in  Sections  3  and  U,  which  are  based  on  the 
representation  (U),  yield 

^  ^        2         8in(2ka-  J)^ 
(15)  (T  ^   n<'.'  (i  -  —      TTT-    )   y  ka  »  1 

as  the  leading  tenos  of  an  asymptotic  development  in  reciprocal  powers  of  ka.     In 
deii.ving  the  frequency-dependent  term  of  (l5),  which  is  believed  to  be   accurate, 
the  normal  derivative  of  the  wave  function  has  the  proper  geometrical  behavior  at 


This  result  was  contained  in  an  abstract  by  the  author  for  the  Symposium  on 
Microwave  Optics,  McGill  University,  Montreal,  Canada,   June  1953. 
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the screen  near  the  rim  (where  it  tends  to  infinity  as  the  inverse  square  root  of 
the  distance),  and  in  the  aperture.  The  importance  of  suitable  behavior  for  the 
screen  distribution,  which  covers  an  infinite  domain,  is  pointed  up  by  inadequacy 
of  the  following  (Kirchhoff  type)  approximation:  The  normal  derivative  at  the 
screen  is  computed  from  the  representation  (3),  when  aperture  and  incident  values 
of  the  wave  function  are  identified.  The  resulting  screen  distribution,  with  a 
stronger  rim  singularity,  will  be  shown  in  Section  U  to  imply  a  deviation  of  0((ka)~  ) 
from  the  geometrical  optics  cross-section.  But  the  geometrical  behavior  near  the 
rim  is  not  sufficient  in  itself  to  ensure  the  proper  deviation  from  the  geometrical 
optics  cross-section,  as  evidenced  by  previously  described  application  of  the  Braun- 
bek  distribution  and  cross-section  theorem. 

Aperture  distributions  at  high  frequencies  can  be  approximated  in  several  ways, 
although  these  do  not  reproduce  the  expression  (l5).  It  is  significant  that  the 
form  of  (15)  can  be  achieved  with  an  aperture  wave  function  that  does  not  vanish 
at  the  rim,  as  is  characteristic  of  the  proper  fvmction  which  varies  locally  as 
the  square  root  of  the  intervening  distance.  This  feature  can  be  verified  in  the 
framework  of  another  Kirchhoff  type  calculation  which  involves  both  the  representa- 
tion (3)  and  (U).  Thus,  from  (3),  it  follows  rigorously  that 


(16)  A  -  -  1^  j  T(p)  dS, 


and  if  now  the  wave  function  (U)  is  approximated  in  the  Kirchhoff  sense  by  confining 
the  integral  to  A,  viz: 


(17)  Tj^(?)  -  -2ik 


) 


( 

G(?,p«)  ds; 


A 


the  result  of  combining  (16)  and  (17)  yields 
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k2      r      ^^ 
(15)  Aj^  .  -  L.         G(p,p')  dSdS'. 


When  the  amplitude  A^     (where  the  index  K  denotes  Kirchhoff  approximation)  is 
employed  in  the  cress-section  theorem,  it  follows  that 

(19)  C^  -  2k  Ira   I     G(p,??')dSdS'. 

^A 

For  a  circvilar  aperture,  the  integral  in  (l9)  is  readily  evaluated,  with  the 
result 

(20)  <r^     -     n,^   (l  .  ^t__J) 


^  y-          ^       sin(2ka  -  j-) 
(21)  iina^fl-i.    tt^   )  »         ka  »  1, 


r        1      sxnuKa--;^ 


Vn         (ka) 


while  the  aperture  wave  function  can  be  expressed  as 


CO 

(22)  T„(p)     -     ka     I  '^         J„(^p)J,(^a)   . 


^  yi:^ 


Thus  our  result  (l5)   for  the  deviation  from  the  geometrical  cross-section 
has  a  coefficient  which  is  twice  that  of  the  foregoing  Kirchhoff  approximation} 
the  associated  wave  function  T„(a)  is  different  from  zero,  and  in  fact  has  the 
value  1/2  if  the  radical  is  replaced  by  k,  this  being  the  mean  of  the  values  1 


*  The  same  aperture  distribution  can  be  obtained  if  the   screen  is  regarded  as  an 
obstacle  in  otherwise  free  space,  with  the  complete  wave  function  a  superposi- 
tion of  incic3ent  and  scattered  parts}  the  discontinuity  in  normal  deirivative  of 
the  wave  function  at  the  screen,  or   'source',   of  the  scattered  field  need  only 
be  approximated  by  its  value  in  the  absence  of  an  aperture. 
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and  0  for  p  <  a,  p  >  a,  respectively.     The  improvement  in  the  cross-section  afforded 
by  the  factor  2,  even  at  moderately  small  values  of  ka,  is  evident  if  we  compare  the 
numerical  values  from  (2i)  vd.th  the  exact  calculations'-  *  J   (see  Table  II).     A  reduc- 
tion of  the  gap  between  (15)   and   (21)  is  achieved  if  we  augment  (l?)  by  a  screen 
integral  with  the  primary  straight-edge  distribution;  it  turns  out  that  the  coeffi- 
cient of  the  second  term  in  (2i)   is  then  increased  by  a  factor  </?  • 

Frahn^    J   has  adopted  Braunbek's  proposal  for  the  aperture  distribution  (which 
vanishes  at  the  irim),  and,  using  (l6)  and  the  cross-section  theorem,  claims  that 


^         ^  sin(ka  +  J)n 

^  *  m  (ka)^/^     ) 


C   2^   na      (  1  +  — i —       ^nrr—  1 »  ka  »  i. 


Vn  (ka)- 

Unfortunately,   there  appears  to  be  an  error  in  the  reduction  of  the  cross-section  to 
its  asyBiptotic  dependence  on  kaj  this  error  invalidates  the  second  term  above,  which 
is  actually  0((ka)      ),  without  the  trigonometric  factor.     Indeed,  for  an  aperture 
distribution  of  this  type,   some  other  means  of  computing  the  cross-section  is  nec- 
essary to  obtain  a  functional  behavior  as  in  (l5)j  one  possibility  is  furnished 
by  a  variational  principle,   as  will  be  demonstrated  in  Section  5. 

Without  entering  into  detail  here,  we  may  appreciate  the  reason  why  high- 
frequency  transmission  characteristics  of  apertures  differ  in  soft  and  hard  screens, 
as  typified  in  Fig.  l.     The  distinction  arises  from  the  differences  in  the  strength 
of  interaction   'across'   the  aperture:  none  at  all  is  provided  by  a  plane  wave   at 
grazing  incidence   (or  even  a  remote  coplanar  line  source)  on  a  semi-infinite  hard 
screen,  whereas  the  same  wave  incident  on  a  soft  sci^en  does  imply  such  interaction. 
Thus,  for  an  aperture  in  a  hard  screen,  the  correction  to  the  geometrical  cross- 
section  represents  a  higher-order  effect  and  appears  with  a  larger  inverse  power 
of  ka  than  if  the  screen  is  soft. 
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Table  II 

IVansmission  coefficient  for  circular  aperture  in  'soft'  screen  for 
normal  incidence* 

t^  '  using  approximate  formula  (i5)« 

(a)  h2l 

t^   using  exact  values  taKen  from  Storruste  and  Bouwkamp'-  ■* . 

ka 

2  06 
2.7 
2.8 
2,9 
3.0 
3.1 
3•■^ 
3.3 
3.U 
3.5 
3.6 
3.7 
3.8 
3.9 
U.O 
U.l 
U.2 
U.3 
U«U 
U.5 
U.6 

Uo7 
U.8 
Uo9 
5.0 
6.0 
7.0 
8,0 
9.0 
10.0 


ka 

,(L) 
*1 

tf) 

.1 

62Uo 

.000003 

.2 

5,7U 

.0000 U86 

.3 

2.27 

.000250 

.U 

.3U9 

.000809 

.5 

.320 

.00203 

.6 

.0220 

.00U36 

.7 

-.m 

.008U3 

.8 

-.1U7 

.0151 

.9 

-.122 

.0255 

1.0 

- .0576 

.O/ill 

1.1 

.0338 

.06u 

1.2 

.lli2 

.097 

1.3 

.261 

.1U2 

1.U 

.385 

.203 

1.5 

.509 

.283 

1.6 

.629 

.38U 

1.7 

.7Uii 

.505 

1.8 

.850 

.6U2 

1.9 

.9U5 

.783 

2.0 

i.029 

.922 

2.1 

i.lOO 

1.0U2 

2.2 

1.158 

1.135 

2.3 

1.202 

1.200 

2.U 

1.233 

1.235 

2.5 

1.251 

1.2U8 

ti- 

tr' 

1.257 

1.2Uii 

1.253 

1.227 

1.2U0 

1.202 

1.218 

1.173 

1.190 

1.1U2 

1.158 

I0I09 

1.122 

1.078 

1.085 

i.0U8 

1.0U8 

1.020 

1.012 

0.99U 

.978 

0.969 

.9U8 

0.9U6 

.923 

.926 

.902 

.908 

.887 

.893 

.877 

.880 

.873 

.872 

.87U 

.867 

.879 

.866 

.889 

.869 

.907 

.878 

.920 

.895 

.939 

.915 

.959 

.937 

.979 

.958 

1.075 

1.065 

.963 

.96U 

.976 

.967 

1.0U2 

1.039 

.987 

.987 
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3*       Screen  integral  equation 

In  this  section  we  base  our  fonnulation  of  the  boundary-value  problem  on  the 
representation  (U).     This  leads  naturally  to  an  integral  equation  for  the  normal 
derivative  of  the  wave  function  at  the  screen*     Designating 

(23)  ik  K(p)     -    -  ^  T(x,y,0),  p  in  S, 

and  recalling  that  3  T(x,y,0)  ■  3  y^  (x,y,0)  in  A,  we  rewrite   (U)   as 


4       Ir  -  p'l  /g  |r  -  p'l 


From  the  requirement  that  the  wave  function  vanish  at  the  screen,  the  basic 
integral  equation  for  K(p)  followsi 

(25)  K(p)  2-3 — 3— dS'-     2__ ;^-dS',       p  in  S. 

-^S        IP  -P'l       4   IP-P'I 

This  equation  holds  for  any  aperture  shape,  but  only  for  normal  plane  wave  excita- 
tion. According  to  (2U),  the  spherical  wave  amplitude  jr   along  the  direction  of 
incidence  is  related  to  K(p)  by 

(26)  y  "^     Ja- J  K(p)  dsj  , 

and  thus  the  aperture  transmission  cross-section  (at  normal  incidence)  becomes 


1" 


(27)  (T  '     A  -  Re  I  K(p)  dS. 

To  improve  on  the  geometrical  value  of  T" ,   a  distribution  K(p)  on  the  screen  is 
needed;  information  about  K(p)  must  evidently  come  from  the  integral  equation* 
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For  an  aperture  of  circular  shape,   and  excitation  by  a  scalar  plane  wave  at 
normal  incidence,   there  is  an  obvious  rotational  symmetry  about  the  normal   (or  z) 
axis  passing  through  the  center  of  the  aperture,  so  that  the  screen  distribution 
is  a  purely  radial  one.     Even  in  this  simple   case   the  integral  equation  cannot  be 
solved  explicitly  by  techniques  now  available   and  we   aim  at  an  approximate   solu- 
tion for  ka  »  1,  i.e.,   aperture  radius  large  compared  to  the  wavelength.     One 
approach  will  now  be  described,  and  another  is  relegated  to  later  consideration 
(Appendix  l). 

Utilizing  the  integral  representation 

iklp-p'l  f°°  yAy  ,     ^^  /7-y      °'       K<i^, 

and  the  addition  formula 

00  •    /     w  \ 

(29)  J,(^|P-P'|)-    H    J„(^P)  J  (^P')e^^^^'^ 

o  *—      lu  in 

-oo 

(where  p,  (jf  are   plane  polar  coordinates  specifying  p),  the  angle  integration  in 
(25)  may  be  performed,  with  the  result 

^00  00  00 


(30)  ^°^       J.(^P)       K(p')J  (^p')p'dp'   "  a  °^         JAKp)JAK^),     P     >     a. 

J      /TT?    0  j  o  ;      ^j-j      0  1 

0  yk  -i^  a  0   yk  -5:^ 


/2 — ? 

Note  that  when  k  ->oo,  and  the  radical  yk  -i^     is  replaced  by  k,  the  left- 
hand  side  of  (30)  is  K(p.)/k,  namely  zero,  since 

(31)  j    K^  JJ^P)JJ^P')     -    ^^-^     . 

o 

Likewise   the  right-hand  side  of  (30)   vanishes  in  this  case,    since 
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.00 

p  >  a. 


00 

(32)  J  d^  J^(^p)J^(^a)  -  0, 


(33) 


More  precisely,  for  ka  »  1,  the  principal  contributions  to  the  ^-integrals  arise 
from  i^  ~  k,  whence  the  Bessel  function  arguments  sire  large  compared  vjith  their 
orders.  Inserting  the  asymptotic  forms 

J^(^p)J^(^p»)  ^  ^  (pp')"^/^(cos  ^(p-p')  +  sin  ^(p+p')'^  , 
J^(^p)J^(j^a)  -  -  ^  (pa)"^/^(sin  ^(p-a)  +  cos  ^(p+a)^  , 


and 


(3U) 


we  obtain  from  (30) 


00  ,00 


(35) 


I  K(p«)[cosJ:(p-p')  +  Bin^(p+p')l  v^  dp' 


.00 


«  _  v/^    ^ sin  ^(p-a)  +  cos  ^(p+a)  ,  p  >  a 

Introducing  the  variable  x  by  a  displacement  transformation 

p  ■  a  +  X, 
and  writing  K(x)  for  K(p),  we  obtain 
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00  00 


yi?7 


0    yk  <        o 
(36)  .00 


I     K(x')   cos  5^(x-x')   +  sin  ^(x+x'+2a)!  >/a+x'     dx' 


-  v^a  — ^       ■         sin  i^  X  +   cos  J^(x+2a)    ,  x  >  0. 


-o  .yi:^ 

Now,   the  integral  on  the  right-hand  side  of  (36)  which  contains  the  sine 
factor  is  perfectly  convergent,  in  fact 

00  kx 

(37)  (        ^— ^    sin  Y^    '     ^     {    H^^^(t)dt,  X  >  0; 


y?I? 


however,   the  integral  with  the  cosine  factor  is  not,  owing  to  the  1/^  singularity 
at  ^  »=  0.     This  defect,   arisirig  from  inadequacy  of  the  asymptotic  forms  (33) >   (3U) 
at  ^  =  0,  is  not  obviated  even  if  we  omit  the  cosine  term  (which  oscillates  more 
rapidly  than  does  the   sine  term  and  may  thus  be  expected  to  yield  a  smaller  rela- 
tive contribution),  since  it  would  then  follow  (via  (37))   that  the  right-hand  side 
of  the  modified  integral  equation  has  the  limiting  value ^  for  k  ->oo,  whence 

K(x)  is  frequency  independent  in  the  limit,   and  8  Y  tends  to  infinity;   this  con- 

z 

flicts  with  the  shadow  notion. 

A  way  out  of  the  difficulty  is  provided  by  first  rewriting  the  inhcmogeneous 
term  of  the  integral  equation  (3)  in  the  numerically  equivalent  form  (see  Eq.  (32)) 


r 

a    AK 
o 


J^iKP)  J^iKa), 


and  then  inserting  asymptotic  expressions  for  the  Bessel  functions.  This  yields 
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(38) 


00  00 

(      ^ 


)/^     '0 


j     I(x')  COS  kU-x^)  +  sin  ^(x+x'+2a)| 


dx« 


oo 


Lyi?:? 


1 1 

~  k 


jsin  Ot  +  cos  ^(x+2a)    ,  x  >  0, 


where 


(39) 


I(x)     =     K(x)    /lT|     . 


It  is  useful  to  note  that 


(UO) 


00        ,- 


-^'~'\ 


.00 


sin  J^x 


di:      e 


iXx 


^   I    ^^rn 


-co 


yi:^ 


6   ->   0,    X   >   Oj 


hence  we  may  recast  (38)   in  the  fonn 


00  .00 

f     <^ 


)/^  i 


j  I(x«)  COS  ^(x-x')   +  sin  ^(x+x«+2a)|   dx' 


(ai) 


.00 


1 
2 


d^        e 


ii^x 


.OO  f- 


-00 


'*''   y^ 


<        'o 


L/1?:? 


1 

k 


cos  J^(x+2a),  X  >  0, 


The  inhomogeneous  terms  are  now  well  defined  and,  as  will  become  clear,  they  have 
suitable  behavior  for  k  — >oo. 

Equation  (hi)  is  fundamental  to  our  calculation  of  transmissLon  through  a 
circular  aperture  in  a  soft  screen.     While  the  equation  does  not  admit  a  solution 
in  closed  form,  it  is  suitable  for  the  purpose  of  establishing  the  first  frequency 
dependent  correction  to  the  geometrical  cross-section,  since  without  the  terms 
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containing  the  rapidly  varying  triponometric  factors     and  the   aperture  raciius  a, 
an  explicit  solution  is  Dossible ,     V/e  recognize  the  solution  as  pertaining  to  the 
shadow  distribution  on  a  half-plane,  located  on  the  positive  x-axis.     With  this 
information  a  perturbation  procedure  is  readily  applied  to   (Ul)   for  the  next  approx- 
imation; this  leads  to  another  half-plane  problem  with  grazing  incidence,   as  dis- 
cussed previously. 

We  consider  first  the   straight  edge  problem  with  normally  incident  excitation, 
using  the  technique  of  complex  Fourier  integrals  to  solve  the  relevant  integral 
equation 

00  00  00  00 

f  _d^_  r  i(o)(^,),,,^(,.,,)d,.  1  r  _^^    i(^'\x.)e^(^-^'^dx. 

(U2) 


f|V^x.)H^'>(k|x-x.|)dx.   -ij      ^^.     ->0. 

O  -00 


y<^ 


Note  that  when  x  >  0, 

!  ^-.ie      ^^-^  ikj      t       ^ ikj     "o     ^^^ 

<.co  i^<  1        \l\r  -1  kx 

(U3) 


,[i(kx  -  J)], 


-    \  /M    «^P|i(kx  -  r)l,         kx  ^  00,     -c 


whereas  for  x  <  C, 

00  .y  00  .j^,  oo 

retaining  e  ><  0  in  the  exponential  in  (UU)  ensures  that  the  entire  function  vanishes 
for  X  -o-oo. 


-  20  - 
Let  us  extend  (U2)  by  defining  an  equation  for  -oo  <  x  <  oo,  namely 


00  00  oo  .^ 

L  yC^  L  loo ' "  •!?:? 


with  the  stipulations 

(U6)  I^°^(x)   =  0,       X  <  0,  F(x)   =  0,  X  >  0, 


It  is  clear  that  for  small  negative  values  of  x,  the  function  F(x)  is  proportional 
to  the  wave  function  off  the  half -plane,  i.e.,  the  aperture  field j  in  fact, 


F(x)  -    /l  +  J     Y(x). 


Multiplying  each  term  in  (U5)  by  e      and  integrating  over  all  x,  we  arrive  at 
the  Fourier  transform  relation  (again  expressed  in  terms  of  a  variable  ^) 


/?T7     ^'"  ^^ 


.00 


where 

(U8)  l(°\^)     .      f    e-^^  I^°\x)dx, 


and 

(a9)  F(^) 


/° 


-00 


e"^^  F(x)dx. 


From  the  obvious  physical  inference  that  I^  '(x)  has  a  propagation  or  wave 
vector  of  magnitude  k  in  the  positive  x-direction,  it  follows  via  (U8)  that  the 
transform  I  '(^)  is  regular  if  Im  ^  <  Im  k  («  e,  say).  The  region  of  regularity 
of  F(^),  determined  by  (Ml)  and  (U5),  in  conjunction  with  the  asymptotic  behavior 
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of  the  Hankel  function  H       (k|x-x'|)  as  x  ->-oo,   is  Im  ^  >  -  Im  k,     V/e  may  thus 
regard  (U7)   as  a  meaningful  relation  in  the   strip    |lin  JC|  <  Im  k  of  the  ^-plane. 
The  transform  equation  (U7)   is  next  rewritten  (in  accordance  with  the  Wiener- 
Kopf  technique)   so  as  to  isolate  terms  regiilar  in  the  upper  (Im  ^  >  -  Im  k)   and 
lower  (Im  ^  <  Im  k)  overlapping  half-planes,   viz: 


(50) 


^+i6 


_/^_. 


/k+ie 


^+ie 


==  *  ^    ^^"^     ^^^^' 


(lower  half -plane) 


>/k+ie 


(upper  half -plane) 


An  entire  function  is  thereby  defined  throughout  the  ^-plane,   and  this  function 
proves  to  be   zero,  by  exaadnation  of  its  asymptotic  behavior  in  the  respective 
half -planes,  using  the  relations  I^°^(^)  -^  (i^)"-'"''^,    \Y,\  ^oo,  Im  ^  <  0,   and 
F(^)  r^  (-ii^)     ,    1^1  ->oo,  Im  J^  >  0.     Thus,  proceeding  to  the  limit,  Im  k  -  e  =  0, 
we  deduce  that 


(51) 


l^°\^) 


l(^-A4) 


(52) 


k 


is  the  transform  solution  of  the  integral  equation  (U2).     By  inversion. 


(53) 


oo 
>-oo 


where  the  path  of  integration  passes  under  the  branch  point  at  5^  -  k.  On  deforming 
the  contour  so  as  to  encircle  a  branch  cut  rxinning  horizontally  to  the  right  from 
^  -  k,  it  follows  readily  that 
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from  which  the  behavior  of  I^  '(x)  for  kx  «  1  and  for  kx  »  1  can  be  inferred. 
The  singularity  at  the  rim  (x  -  0)  is  of  a  familiar  type. 

Using  this  solution,  we  may  inquire  about  the  corresponding  scattered  ampli- 
tude ^         ,   given  in  Eq.  (26),  which  we  now  cast  in  the  form 


(55)       j{<-°'>    .-g 


na  -  2na 


p'\.)/r:i  dxl 

o  •i 


If  the  radical  is  approximated  by  1  +  •=—  ,  we  encounter  the  integrals  (cf»(52)) 

(56)         ri(°^x)dx  .  l(°)(0)  -  1^ 


and 


,CD 


8k' 


(57)  xl(°)(x)dx  =  il(0)'(0).--^- 


whence 


f 


(58)  ;<(0).   -i^-l-j^,  k.»l. 


A  comparison  with  Eq.   (12)  reveals  that  the  term  -a/2  is  indeed  the  Braunbek 
correction}  the  i/k  term  gives  a  contribution  of  order  1/k     to  the  cross-section, 

but  as  we  will  show  presently,  a  better  approximation  to  the  solution  of  the  inte- 

-3/2 
gral  equation  (W.)  generates  a  lower-order  contribution,  namely  k     '     in  C^.     The 
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ease  vdth  which  the   scattered  amplitude  follows,   once   the   transfcnn  of  the   screen 
distribution  is  known,  will  be  noted. 

The  Fourier  integral  technique  is  not  as  successful  in  coping  with  the  inte- 
gral equation  (hi)  for  l(x),   since  the  presence  of  a  terra  involving  x  4-  x'   in  the 
kernel  gives  rise  to  an  integral  equation,  rather  than  an  algebraic  relation,  in 
the  transform  variable.     The  latter  integral  equation  can  be  solved  by  iteration; 
but  we  shall  apply  instead  a  direct  perturbation  procedure  to  (Ul),  which  will  give 
the  desired  information  equally  well.     Let  us  therefore  write 


(59) 


I(x)     -     l^^hx)     +     I^^^(x), 


where  I       (x)  is  the  function  defined  by  (Ii2),  whose  transform  is  contained  in 
(51),  while  I     ^(x)  is  defined  by  the  integral  equation 


-CO 


.00 


.00  I- 


^ 


^-      l(l>(x.)e^^^-^'^dx. 


T-l 


-00  yk  -J^     o 


Ly(?7 


COS  l^(x+2a) 


(60) 
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I^    ''(x«)sin  i:(x+x'+2a)dx',     x  >  0, 
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The  first  integral  on  the  right-hand  side  of  (60)  can  be  transformed  according  to 
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Lyk^ 
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cos  ^(x+2a) 


?  S^(t)   -  i  J  J,(t)   -  i   j.    dt 


:(x+2a) 


(61) 
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^  V  I         Z~~Z ^•'-  '  ■§  ^  k(x+2a) 

/k(x+2a) 


.  +  Onc(x+2a)  V   ,  ka  »  1, 
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where  S^  denotes  the  Struve  function  of  zero  order.   Instead  of  discussing  the 
second  integral  in  (60)  separately,  it  is  convenient  to  combine  the  two,  whence, 
with  the  help  of  (1^8)  and  (5l),  we  find 
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'-00 


d^ 


v^:v  'o 


[  Ji)f   ,.  i^(x-x'),  , 


i^(x+2a) 


i""^)  .  1 


V^ 


'V^^ 


-i^(x+2a) 


(62) 


1^°^^)  ±1      1 


1  r^   f  i.(x.2a)(^  .  ^^  ^   ^.i.(x.2a)(^  .  ^0  ^  ^  ,  ,, 


This  integral  equation  -  which  differs  from  (U2)  in  the  structure  of  the 

inhomogeneous  term  -  also  describes  a  half-plane  diffraction  problem.  Since 

I^   '(x)  is  viewed  as  a  perturbation  on  I^     (x),  we  may  avail  ourselves  of  a 

simplification,   appropriate  for  ka  »  1,  in  the  forcing  terms  of  (62),  Now, 
for  a  »  1 


(63a) 


[>»^-(-;^)-^/Fe4u-5,].^.0(^) 


la        * 
> 


Ihe  leading  term  of  (63a)  is  easily  obtained  by  recognizing  that  the  integral  de- 
pends  principally  on  values  of  tSJl  when  a  »   1.  Hence,  writing  1-t  =  u,  discarding 
the  constant  relative  to  the  radical,  and  placing  the  reciprocal  of  t  equal  to  unity, 
we  find  (see  (93)) 

"00  "00 
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(63b) 


dt     ^-iat 


(?":^y  ^ 
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SO  that  on  retaining  only  the  loirest-order  term  for  (62),  we  have 


(6U) 
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iooVCV^ 
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It  is  apparent  that  (6U)  refers  to  a  half -plane  excited  by  a  remote  line 
source  which  is  parallel  to  and  coplanar  with  its  edge.     In  the  present  context, 
this  source  localizes  the  total  effect  of  the  far  side  of  the  screen  and  the 
aperture  as  regards  interaction  in  the  plane  of  the   screen.     For  the  purpose  of 
asceirtaining  the  function  I       (x)   to  the  lowest  order  in  1/ka,  the  line-source 
excitation  in  {bh)  can  be  replaced  by  a  plane-wave  excitation,  namely 


,00 
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.00 

j     l(-^\x.)e^^^"-"'^dx.   =  -  l/^  exp[iic(x.2a)-  ^],  x  >  0. 


The  integral  equations  (6U)   and  (65)  are   amenable  to  the  transform  technique 
employed  earlier.     For  the   simpler  of  the  two  cases,  involving  (65),  v«  start 
with  an  extended  equation 
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where 


(67) 


I^-'-^x)     -     0,       X  <  0, 


G(x)     =0,       X  >  0. 


Ihe  transform  equation  thus  becomes 
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(68)      _-^^»-j-y__   -__ +  ^  G(^), 
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dx 
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regtilar  for  Ira  ^  >  -  Im  k; 


Im  i^  <  Im  k   , 


Rearranging  Eq.  (68),  defined  in  the  strip  |lm  ^ |  <  Im  k,  so  as  to  separate 
terms  regular  in  the  half -planes  Im  ^  <  Im  k  and  Im  ^  >  -  Im  k,  we  have 
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Since  1^    ' {Y,)  and  G(^)  have  the  same  asymptotic  behavior  as  I^^(^)  and  F(^)   (cf. 
(50)),  the  entire  function  vanishes  and  consequently 
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The  approach  to  Eq«    (6li)   is  similar,   and  instead  of   (66)  we  find 

^^^^(^^  1  r  nn      r     e^(^<^^  1 

(75)  ^,  ^     exp[i(2ka-  i)]  ^  dx+i-G(^). 

Integration  by  parts  shows  that 
'°°  i(k<)x 


(76) 
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rH      ^    ^-^       2(2a)^/2  (kO^    ^     / 


the  first  term  of  this  equation  leads  back  directly  to  (68).  If  the  second  term 
is  also  retained,  the  transform  I^  (X.)   tvirns  out  to  be 


+ 


-00 


where  the  path  of  integration  is  a  line  in  the  strip   |lm  t|  <  Im  k,   and  the  integral 
represents   (proportionally)  that  part  of    /k+Ji  /  (k-^)     which  is  regular  in  the 
half -plane  Im  ^  <  Ini  k«     The  integral  may  be  evaluated  after  con5)letion  of  a  circuit 
in  the  upper  half  of  the  t-plane;  we  have 

tt\  11  r  "li  T  1  i(2ka-  1-)     3   ♦  t 

(78)  I^^^^)  -  ^    -^    exp  i(2ka.  J)     -i- ^  — i^  e  ^  ^, 

^^    v^JKiT  L  ^  J  XT?      16k  A  (ka)^/2  (1-  |)^/2 

and  it  is  verified  that  the  deviation  from  (73)  is  of  higher  order  in  l/ka. 
When  I^   '(x)  is  derived  from  (73)  and  combined  with  I^^(x),   we  find 
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f^^  /,  X     inA  r     ^    i(2ka-  JL)  ikx    C°°     it   -i 

(79)  I(x)  -  l(°>(x)*l(l^x).  ^—  [(1-2  Vg  e      ^  )  2—  ^i    2__  dt]  , 

yiT   I-      ^^  ^    {^  vt  J 

which  shows  a  change  in  the  coefficient  of  the  dominant  term  near  the  rim.  Evalua- 
tion of  the  scattered  amplitude  as  in  (55)- (57),  with  the  inclusion  of  I^  (0)  from 
(7U),  yields 

(80)  A^^^  ^  -  i^  -  I  +  v^  exp[i(2ka-  J)]  +  0(|)  , 


and  thus 

^       sin(2ka-  J)  s 
(81)     <r  -  na'^  (  1  -  i 57§-  )  +  0(ka)  ^  . 

We  have  thus  succeeded  in  deriving  an  explicit  correction  to  the  geometrical 
cross-section,  which  is  of  lower  order  in  lA^^  than  that  obtainable  from  the  pre- 
vious approximation,  and  3hoiJ.d  contain  the  principal  diffraction  effect  when 
ka  »   1. 

U.  Variational  calculation 

It  is  natural  to  consider  the  applicability  of  variational  techniques  for  the 
problem  at  hand,  in  view  of  their  success  for  other  aspects  of  scattering  phenomena. 
Indeed,  the  transmission  through  a  circular  aperture  in  a  soft  screen  has  been  cal- 
culated with  high  accuracy  for  the  low-frequency  range,  0  <  ka  <  U,  by  use  of  a 
variational  formulation  which  involves  approximation  to  the  aperture  distribution 
of  tI-  -I  ♦  Another  variational  principle  was  developed  by  Levine  and  Schwinger'-  ■•, 
^ich  is  closer  to  the  framework  described  in  Section  3  in  that  it  relates  to  the 

screen  distribution  of  3  Tj  this  principle  also  has  the  merit  that  it  incorporates 

z 

the  refined  Kirchhoff  expression  (18)  as  part  of  the  scattering  amplitude.  The 
efficiency  of  variational  procedures  is  contingent  on  the  choice  of  trial  functions. 
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and  here  the  dictates  of  similarity  to  the  proper  distribution  and  relative  simplicity 
for  computational  purposes  nrast  be  considered.  That  these  are  not  always  compatible 
is  evidenced  with  half-plane  screen  distributions,  which  prove  too  complicated  for 
exact  evaluation  of  the  integrals  concerned.  Thus,  in  earlier  application,  only 
simple  forms  of  the  screen  distribution  were  employed,  and  these  preclude  an  accurate 
estimate  of  the  departvtre  from  geometrical  behavior.  It  will  be  our  purpose  here  to 
show  that  more  precise,  though  complicated,  distributions  such  as  that  of  the  half- 
plane  may  be  exploited,  at  the  expense  of  a  consistent  approximation  to  the  integrals. 
We  will,  in  fact,  confimi  by  this  approach  the  explicit  correction  to  the  geometrical 
cross-section  given  by  (8l), 

A  variational  characterization  of  the  scattered  amplitude  A  ,   based  on  the 
integral  equation  (25)>  can  be  extracted  from  the  development  in  (lO],  where  the 
directions  of  excitation  and  observation,  as  well  as  the  shape  of  aperture,  are 
left  arbitrary.  The  result  is 
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with  a  Kirchhoff  term  A^  of  the  form  (18).  The  cross-section  deduced  from  (82), 


(63)  <r-  (Tj^*  ^  Ira 
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likewise  involves  a  Kirchhoff  part,  evaluated  previously  for  the  circular  aperture. 
If  the  stationary  nature  of  (83)  is  utilized  in  effecting  a  scale  transformation 
of  K(p),  we  obtain  the  cross-section  in  a  homogeneous  form: 
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To  improve  the  tractability  of  the  integrals  in  (83)  and  (8U),   we  modify  the 
Green's  function  along  lines  suggested  in  the  preceding  section.     Thus^  f or  a 
circular  aperture,   recalling  (28),    (32),    (33),    (39),   and   (UO),  we  put 
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Similarly, 
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A  notable  feature  of  this  rearrangement  is  that  the  transforms  of  the  screai 
distribution  are  put  in  evidence,  as  these  take  a  simpler  form  in  the  half -plane 
problem  than  the  coordinate  functions.  Using  the  transform  I^  (^),  which  corres- 
ponds to  incident  wave  excitation  (with  scale  fixed  by  (51)),  as  an  initial  approx- 
imation, we  have 
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and 
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Hence,  the  corresponding  cross-section  (83)  becOTies 
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If  the  scale  factor  in  I    were  not  known,  we  should  appeal  to  the  homogeneous 
form  (8U)  of  the  cross-section,  and  deduce  that 
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(92) 
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\Vhen  a  »  1,  the  values  of  I^(a)  and  I- (a)  are  determined  principally  at 
X  :il  in  the  integral  representations  (90),  (91).  Iifriting  x  «  1  -  •■  ,  we  have  for 
a  »1 
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and 


I  (a)  ^   i 


,00 


"00 


^   Jexp[ia(l.  |)J(i.  y?)2  .  exp[-ia(l.  |)J 


>  • 


Since  arg  V^  equals  n/2  or  0  according  as  t  <  0  or  t  >  0,  we  find  for  o  >  0 
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so  that,  for  a  »  1 
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and 

(96)  IgCa)  ^  (  V^  -1)^  /?  exp[i(a-  J)]„ 
From  the  estimates  (95)  and  (96),  we  find 

(97)  2I^(a)  -  l^{a)   ^  -  /^  exp[i(a-  J)j,  a  »  1, 
and  thus  (89)  predicts 

(98)  '^  -  ^K  ■  I  V^  sin(2ka  -  J),  ka  »  1. 

Recalling  the  high-frequency  development  of  (T^  ^Eq. (21 )u  it  is  observed  that  the 
variational  contribution  to  (98)  equals  that  of  the  Kirchhoff  term  of  the  same 
order,  and  so 

/     siji(2ka-  ^)\ 

(99)  C^^na^  1-  4.  ,      ^    L  ka  »  1, 

\  VK       (ka)3/2    J 

which  agrees  with  the  explicit  terras  in  (81).  The  same  asymptotic  behavior  is 
arrived  at  by  employing  the  homogeneous  fonn  (92)  of  the  cross-section,  since  its 
expansion  in  powers  of  L.  and  Ip,  which  functions  approach  zero  with  increasing 
argument,  coimencos  with  the  linear  combination  21,  -  !„  appearing  in  (89) • 

A  detailed  investigation  of  the  behavior  of  I^  and  Ip  for  large  and 
small  values  of  the  argument  is  feasible,  although  we  can  circumvent  this  consider- 
ably in  our  fvurther  remarks  about  the  related  cross-sections.  Indeed,  since 
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the  combination  21,  -  I„  appearing  in  the  cross-section  (89)  is  related  to  a  new 
pair  of  integrals,  one  of  which  has  previously  been  studied  (cf.  (6l)).  The  lead- 
ing term  in  the  asymptotic  development  of  this  integral  (a  »  1)  suffices  for 
reproducing  (97),  and  it  is  easy  to  verify  that  the  residual  term  in  the  cross- 
section  (99),  stemming  from  the  other  integral,  is  Ouka)"  j.  The  relation  (100) 
is  also  convenient  for  ascertaining  the  low-frequency  (ka  «  l)  behavior  of  the 
cross-section;  we  note  simply  that  whereas  the  transmission  coefficient  (cross-section 
relative  to  apertiu?e  area)  should  follow  the  Rayleigh  law  and  decrease  as  (ka)  ,  the 
forms  (89)  and  (92)  tend  to  infinity  as  ka  ♦  0. 

Ue  next  examine  the  effect  of  including  a  secondary  (grazing)  wave  terra 
in  the  transform  I(^),  If  I(^)  -  I^°^(^)  +  l'"'"^(^),  then  possible  additional  correc- 
tions to  the  geometrical  cross-section  of  the  order  provided  by  I^   alone  (cf.  (99)) 
are  to  be  sought  in  the  leading  terms  of  (85)  and  (86)   We  find 
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where  the  paths  of  integration  are  lines  in  the  strip  |lm  ^|  <  Im  k.     Utilizing  the 
explicit  forms  of  I^     ,    I^^    (the  latter  to  within  a  scale  factor),  we  obtain 
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where  the  crucial  stage  is  reached  after  deformation  of  the  integration  contours  in 
the  upper  half  of  the  ^-plane.  Thus,  I^   does  not  contribute  to  the  explicit  correc- 
tion term  of  (99)  in  the  variational  calculation,  although  as  we  have  seen  previously, 
it  yields  the  entire  contribution  in  a  non-variational  procedure.  The  variational 
technique  affords  us  special  advantage  here,  as  the  integral  equation  may  be  solved 
to  the  requisite  order  of  approximation;  however,  the  latter  situation  does  not  occur 
with  the  aperture  formulation  given  in  the  following  section,  and  then  a  varia- 
tional procedure  is  appropriate. 

Before  concluding  with  screen  formulations, we  shall  consider  in  more  detail 
the  KLrchhoff  approach,  briefly  alluded  to  in  Section  2,  which  is,  in  a  sense,  recip- 
rocal to  that  discussed  earlier  (cf,  Eq.  (l7)-(22)).  Once  again,  the  representations 
(3)  and  (U)  are  combined,  but  now  (3)  supplies  an  approximate  wave  function,  in  terms 
of  the  ixnperturbed  aperture  distribution,  from  which  a  normal  derivative  at  the  screen 
is  determined;  this  normal  derivative,  together  with  (U),  gives  the  transmitted  field. 
We  write 
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and  (cf.  (23)) 
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using  the  differential  equation 
(103) 
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and  the  relation 
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With  use  of  the  integral  representation  (28),  it  follows  that  in  the  case  of  a 
circular  aperture. 


(105) 


where 


«p)  '-li^^hh^  ^'y^. 


.00 


F(p) 


d^ 


^/^ 


J^(^p)J^(^a). 


For  large  ka  we  find  (see  Section  3)  that  the  function  F(p)  is  given  approximately  by 

,00 

(106) 


F(P)^  k    ^    i  H^^\t)dt, 

^  k(p-a) 

and  the  dominant  singularity  of  K(p),  p  ♦  a,  is  therefore  of  the  form  (p-a)  ", 
with  imaginary  coefficient.  Since  only  the  real  part  of  K(p)  contributes  to  the 
cross-section. 


v-T 
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(107)         <r-  na^  -  2n  Re  I  K(p)pdp, 


/" 


an  infinite  value  of  <r  is  avoided.  For  the  evaluation  of  (107 ),  the  real  part 

of  F  is  significant,  namely 

k 
Re  F(p)  -  f    /^    J  (^p)J  (^a), 

and  it  may  be  shown  that 

2ka  J  ,s 

(108)  Cr-  na^  *^    j     -LL  dx  ^  na^d  ♦  ^),      k  »  1. 

'o 

The  excessive  magnitude  of  this  correction  to  the  geometrical  cross-section,  for 
ka  »  1,  also  occurs  in  an  approach  to  be  described  at  the  end  of  Section  5» 

5.  Aperture  integral  equation 

From  the  expression  (16)  for  the  scattered  amplitude  ^  and  the  general  theorem 
(13),  we  obtain 

(109)  <y  '     Re  J  I(p)  dS 

as  the  cross-section  of  an  aperture  in  a  soft  screen  at  normal  incidence.   If  we 
put 

(no)      Y(?)  -  1  +  -/C  (p), 

the  geometrical  part  of  the  cross-section  is  isolated,  and  we  have 

(111)         (T-  A  +  Re  J  'X(?)dS, 

A 

the  analogue  of  (27). 
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One  might  think  that  the  aperture  distribution  f (p),  defined  in  a  bounded 
region,  would  be  more  easily  characterized  than  a  screen  distribution  K(p),  which 
covers  an  infinite  region,  and  indeed  for  low  frequencies  this  is  the  case.  How- 
ever, in  the  high-frequency  case,  the  aperture  field  determination  is  more  diffi- 
cult, particularly  when  it  is  desired  to  account  for  the  behavior  both  near  the 
rim  and  near  the  aperture  center.  The  field  distribution  near  the  aperture  rim 
is  readily  obtained  for  the  simple  strain; ht. edge  geometry,  in  which  case  the  screen 
and  aperture  problems  are  complementaiy.  For  the  circular  aperture,  however,  the 
use  of  a  straight-edge  approximation  is  necessarily  inaccurate  near  the  center; 
for  better  approximations,  the  complementarity  is  entirely  lost«  It  is  somewhat 
fortiinate  for  this  approach  that  the  parts  of  the  aperture  distributions  ¥  and 
^]f\,   near  the  rim  are  weighted  more  heavily  in  the  cross-section  calculation. 

To  bear  cut  these  contentions,  and  also  to  find  the  quantitative  predictions 
of  aperture  formulations,  we  commence  with  the  exact  integro-diff erential  equation 
in  the  aperture,  which  is  a  consequence  of  (3)  and  the  property  that  in  A 

z       z        • 
This  integro-differential  equation  is 

9    0     f      ^        oik|p-p'| 
(112)         (7^  +  k*^)  I  T(p')  2-_ — -—dS-  -  2iiik,  p  in  A  . 

/.        I P  ■•  P '  I 


A 


Here  V"  is  the  transverse  (to  z)  part  of  the  Laplacian  operator.  When  the  aperture 


t 
is  circular,  (ll2)  can  be  written 


2  r°°  (^ 

Z      p  3p    •  I    rfK — T  o' 
o 


(113)       (^*i  h'^^A     -7£=:V^P^ 


p'J^(^P')T(p')dp'  .  k,    p  <  a 
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and  integration  yields  an  ordinary  integral  equation 


.00 


(liU) 


d^ 


/o     \/k^-^ 


J^iKp)   I   p'JQ(^p')T(p')dp'  -  i  +  aj^(kp)   ♦  pN^(kp),        p  <  a. 


From  the  condition  of  finiteness  at  p  ■  0,  we  may  require  P  ■  0;  the  remaining  in- 
tegration constant  a  is  chosen  to  ensure  that  1(a)  =  0.  The  function  "^  (p) feq. 
(110))  is  described  by  the  equation 


.00 


(115) 


vC^    ° 


r 


^^—    JJKp)   I  P'J„(^P')  X(p')dp' 


,00 


aj^(kp)  *  a  I  dS  ^i  -   1 


y^^ 


•'ofep'^'ltea). 


since 


,00 


(116)         j  Jq(^p)  J3^(^a)  d^  -  I  , 
'o 

and  7^ (p)  is  subject  to  the  boundary  value 


p  <  a. 


p  <  a 


(117) 


X(a)  -  -1. 


Placing  p  ■  0  in  (ll5)  relates  a  with  an  integral  of    7((p);  we  have 


(118) 


a  -  -a  1     d^  Ji  -  ~f~  }  J  (^a)  ♦ 


00  a 


yi?^ 


pj  (j;p)X(p)dp 


^ika.^         e^^PX(p)dp, 
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and  this  can  be  used  as  a  consistency  check  on  the  value  of  ot  obtained  via  the 
boundary  condition  (117 )• 

V/hen  ka  »  1,  the  important  range  is  ^  :^  k  in  the  integrals  of  (ll5),  and 
the  Bessel  function  arguments  are  therefore  large  in  comparison  with  the  orders, 
except  near  the  aperture  center  (p  ■  0),  Disregarding  the  latter  aspect,  and 
approximating  to  the  Bessel  functions  by  asymptotic  forms  in  the  familiar  manner, 
we  find 


00 

d^ 


X(p')  V^  [cos  J^(p-p')+sin  ^(p+p')]dp'  "  na  y^  cos(kp  -  J) 
^o  Vk*-<~  'o 

(119) 


^/?:? 


V^   I  X   [  ,^     m    W^  ^(p-a)+cos  ^(p+a)j, 


o    Vk-C  o 

00 

+  -/a   I   -^1   ~    -  ^1  |sin  ^^p-a;+cos  K{p*a)\,  p  <  a. 

'o      /k^<" 


Introduce  the  coordinates  x,x*  such  that  p"a-x,  p'-a-x',  and  define 

(120)  J^(X)  -  y/i  -  I     X(X), 

whence 


^00         a 

d^ 


4  yjT-: 


I  0(x')rcos  ^(x-x»)+sin  ^(2a-x-x»)|  <ix« 


K      '0 
(121)  00 

"  ''/h  cos(k(a-x)-  £■)  "  I  T  [— ^'— •  -  ^[^^^  -  cos^(2a-x)],   X  >  0. 

This  integral  equation  holds  when  x  is  positive  and  less  than  a.  If  we  omit 
the  integrands  which  contain  the  aperture  radius  a,  and  take  the  upper  limit  of 
the  x'  integral  as  infinite,  the  corresponding  function  5^(x),  0  <  x  <  oo,  can  be 


-  Ul  - 


obtained  exactly,  ^(x)  will  deviate  from  j^(x)  near  x  ■  a,  the  aperture  center, 
and  actually  (121)  itself  is  not  appropriate  in  this  region,  but  we  may  still 
expect  to  find  a  reasonable  cross-section,  on  the  ground  that  the  part  of  the 
aperture  field  near  the  rim  (x»0)  is  more  important. 

In  the  approximation  contemplated,  the  functicn  ^(x)  has  a  local  character, 
and  consequently  on  resolving  the  trigonometric  term  (with  a  as  a  coefficient) 
into  plane  waves,  it  is  clear  that  the  wave  e    ,  describing  interactions  from 
distant  parts  of  the  aperture,  must  be  without  effect.  We  shall  therefore  omit 
this  forcing  term  (this  is  justified  in  (Appendix  2n.     Consider,  then^the  integral 
equation 


1 


,00 


'-00 


.00 

'^         I       ^f   ,.   ii^(x-x«),  . 
J^     I   $2^(x')e       dx- 

•^     '-oo 


""   \/m  ^^^   exp[-i(ka-  f)j. 


X  >  0 


,  X  <  0 


.00 


^-00 


d^    e 


i^x 


>/^ 


+  /^(x),   -CD<X<00 


with  the  stipulations 


(123) 


^(x)  -  0,   x  <  Oj    A  (x)  -  0,  X  >  0. 


The  transfora  equation 


(12U) 


J|.a^exp[.i(.a.j],4 


/?:? 


A(^), 


where 


(125) 


.00 


^{K) 


e"^^  pr(x)dx. 
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(126) 


00 

A(^)  -  j  e"^'^  A(x)dx, 
/-oo 


holds  in  the  strip  |lm  ^|  <  Im.  k("6)  of  the  i^-plane.     When  Eq.  (12U)  is  rewritten  as 


2^- 


la 


f2Tika. 


v^ 


-i(ka-  f) 


k< 


i 
:+ie 


_  Vl?<        >/icfi£_ 


(127) 


io 


^2Jika 


exp[-i(>ca-  f,]   .M^ 


the  left- and  right-hand  members  are  regular  in  the  half -planes  Ira^<Imk,   Ira^>-Imk 
respectively,   and  we  deduce  from  the  vanishing  of  the  entire  function  thus  defined 
that 


(128) 


jy(^)-^(i.vC|)>4^  i 


-i(ka-  ?) 


yita        /k  -  J^ 


with  the  inverse 


(129) 


^(x) 


in/U    itac 


+  i 


r°°   it 

2_dtl-    '^ 


^n     /H    > 


exp 


ik(x-a)+  ^ 


kx 


V^   1    n  ya 


•^ 


-1/2  * 

The  terms  involving  x  '  have  to  cancel  and  therefore 


(130) 

whence 

(131) 


-/^ 


i(ka.  I) 


,00 


Ttf   \        i  inA  I   e   J*    1  (   dt 
^(x)  -  —  e  '      dt  -  -  -     -T- 


,00 


ikxt 


^/TT 


Note  that  with  this  value  of  a,  the  corresponding  inhoraogeneous  term  in  the  integral 
equation  (122)  is  independent  of  the  aperture  radius* 
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and 

(132)  ?to)  - 1  (1-  -ir  )  • 

At  X  -  0  (p  ■  a),  both  ^  and  Ji   (obtained  from  (120))  take  the  value  -1,  in 
accordance  with  the  boundary  relation  (117);  for  x  >  0,  ^  is  the  local  straight- 
edge function,  decreasing  as  kx  becomes  larger,  and  due  to  the  factor  yl-  —  which 
links  5^  and  yL  ,   the  latter  tends  to  infinity  at  the  center  of  the  aperture. 

The  expression  (132)  for  ^(K)   can  in  fact  be  derived  from  the  transform  solu- 
tion of  the  integral  equation  (U5)>  since,  as  already  pointed  out,  the  function 
F(x)  appearing  therein  is  related  to  T(x),  namely 

(133)         F(x)  -  -  ^   y/fTf  Kx),      X  <  0. 


Now 


2jlL  __1 

k^ 


F(^)  -  -  ,,   , 


^ 


SO  that  the  tranaf onn  of  yl+  —  T(x)  is  rr  — ——--  »  and  subtracting  the  transform 
h  of  the  incident  field  in  -co  <  x  <  0,  we  find 

^(^)  -  r  i-^ —  - 1); 


K 


v^l 


this  is  substantially  the  result  (132),  except  for  reversal  in  sign  of  ^,  which 

stems  from  the  fact  that  here  the  aperture  domain  corresponds  to  negative  x. 

The  aj>5)litude  ^^  can  be  written 

a 


-  uu  - 

and  if  the  radical  is  approximated  ^7  1  ~  TrT  »  ^^'^   'the  integration  limit  there- 
after is  set  equal  to  infinity,  then  usiny  (132)  we  find  (cf,  (58)) 

(135)        -^"^  -^-?*lik'       ka»l, 

which  implies  a  correction  of  order  (ka)   to  the  geometrical  cross-section.  The 

retention  of  a  finite  upper  limit  in  this  integration  has  the  effect  of  replacing 

—2  —3  /2 

the  (ka)   correction  with  one  of  order  (ka)  '  ,  where,  however,  the  accompanying 

trigonometric  factor  depends  on  ka  and  not  2ka,  Actually,  the  radical  in  (13U) 

need  not  be  approximated,  and  it  can  be  verified  that  the  cross-section  resulting 

from  the  distribution  ]?(x)  is,  asymptotically. 


^  ^  na^  (1  -  2i2J^),       ka  »  1. 
(ka)2 

It  is  interesting  to  ascertain  the  value  of  o  which  follows  from  the  relation 
(118)  when  the  function  ^(x)  is  involved;  this  value  is 


1   ika  .  ika  |    dx    -ikx 
a  -  -  T-  e   -1  e     |  e 


3f(x) 

o  /^ 


a 


00 

1   ika  ^  i  ika  (  du    1    f    dx    ik(u-l)x 


(    du    1 


e    +  —  e    I  -r I  137  e 

a 


For  X  j(  0   and  kx  »  1,  asymptotic  evaluation  of  the  u  integral  yields  yrrr  e  '  , 


whence 


a    -•/Tf    expfldta- J)j, 


a' 


ka  »  1 

K        -  L  u j- 

which  agrees  with  (130). 
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Further  progress  with  the  integral  equations  for  V    or  <if  Is   difficult;  an 
estimate  has  been  made  of  the  change  brought  about  by  extending  the  integration 
limit,  whereby  J^  ♦  ]?,  but  this  did  not  significantly  improve  the  cross-section. 
We  shall  therefore  turn  next  to  a  variational  formulaticai  of  the  scattered  ampli- 
tude, and  see  how  our  present  aperture  distribution  can  be  exploited.  The  varia- 
tional expression  reads 

^  -  -  ij^    T(?)dS  ♦  (^)    I(p)  2__-__T(p.)dSdS' 
A  /i       I P-P '  I 


(136) 


(2n) 


i-^  I  nit)   '^'fCp') 


,ik|p-p' 
Ip-P'l 


dSdS', 


with  aperture  functions  Y  that  are  required  to  vanish  at  the  rim.  Writing  T  ■  1+X , 
we  have 


V--^(A+  )/(?)dS)+(|^) 


Ik  I P-P 'I 


ik|p-p'i 


® dSdS-  +  2   ^  '^^         X(p')dSd3t 


IP-P' 


'A 


|p-p' 


(137) 


A       I P-P' I 


--1-^   7/(?).7«  (?.)  ^ 


ik|p-p'| 


(2n) 


I P-P' I 


whence 


(T-  2A  +  2Re  I  /(p)dS-  ^   Im 


f  ik|p-p'|        /  ik|p-p'| 
^  _^_^ dSdS«+  2  2__-__ ;i(pi)dSdS« 

l4    Ip-p'I  a    I'^p'I 


(138) 


f      ik|p-p', 
J/(p)  V^..  ■   ^(p')dSdS' 


|p-p'l 


"I     1  ik|p^p'| 

>  ^  Im   v/(p).V.X(p.)  2___^dSd3' 


2nk 


P-P'l 
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We  now  particularize  the  various  integrals  in  (138)  containing  the  Green's  func- 
tion, for  a  circular  domain  A  and  a  purely  radial  function  ^(p).  The  first  may 
be  handled  exactly  (cf.  (19),  (20)): 


(139) 


k  T 


(    ^ik|p-p«|  JT(2ka) 

-^ dSdS'  -  -Tta^d-  ^  ^     ) 


A   I  p-p 


-o> 


^   sin(2ka-  ^) 
o<  -na(l rjo —  » 


ka  »  1; 


and  in  the  other  integrals  the  familiar  Bessel  function  approximations  are  utilized. 
Thus 


ik I p-p' 


^'^■^•1 V-  2r"      1      1 

^^-^ X(p')dSdS'  -  ia(2n)'^       d;^  ^        "^         -  ^V  J,  (^a) 


j  pX(p)J^(^p}dp  +i^^    J   p?^(p)dp 
o  ^ 


,00 


00 


-«a  ^    2iSi_  -  ,„.a  (    f 


-00 


v<?:? 


1       1 


[^'  i 


2i^a, 


e^^V(^)  *  e-'^^"j2((-^) 


-2i^a, 


+  il2^    [p/(p)dp. 


where 


(lUl) 


0{K)  -       e'^^  0(x)dx 


Note  that  the  integral  of  ^  over  the  aperture  is  such  as  to  cancel  with  the 
second  term  in  (138).  Further, 
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JX(P)    ^^  r^       ;C(?r)dSdS'   -   i(2n)2   j      -^  [  p/Cp)  J^(^p)dp    fp '/(p- )  J^(^p' )dp  ' 

4  |p-p'l  A  VkW^  ^ 


(lli2) 


ii  2id.a 


^00  00 

^ tif'^\rif.^\^r,„^     I       ^ 


0(^)0«)+2na 


^cx.    yi?<^ 


v^ 


2i^a, 


(0(^2.  ^-2iJ^a  ^^_^^y 


The  remaining  integral,   before  approximatior^  is 


j  vX(p)  •  V'X(p') 


jikl^p'l 


dSdS» 


00  (^  f 

(1U3)       -  i  E  (2-6)„„)   I     -^^     I  dSdS-  ^t2L  cos(0-jZi')J  (^p)J  (^p')cos  n(0-{2f«) 


-  i(2n) 


,00  .a 

2   /  ^d^ 


/?:?/o 


a 


and  to  avoid  a  singularity  in  the  radial  integrals  at  the  lower  limit,   the  insertion 
of  an  asymptotic  form  of  the  Bessel  function  is  deferred  until  after  integration  by- 
parts,   using  the  boundary  value     X  (^)  "  ~1>   ^z: 


a 

^  pJ^(^p)dp  -  -aj^(^a)  -  ^  I    <(p)pj^(^p)dp 

o 

a 
^  -/g  sin(^a  -  J)  -/f     f  X(P)  /?  cos(^p-  J)dp 
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Thus, 


gpj^(^p)dpr  ^n 


exp[2i(^a-  f)](j2^(ii)-  \J*2(^i.X,)'  |)  (0«)*  \) 

+  exp[-2i(^a-  J)](^«K|)^ 


and 


VA(p)   .   V'X(p') dSdS' 


(lUU) 


00 


/-oo  yk"^-^ 


na 


/°°     2 

/o  yi?:? 


.2^^(0(,).  i  )^      .  e-2^^(0(-,).  if  I  . 


Collecting  restilts,   we  have 


..  „2  ^.    1      °^"''^''^"  f' 


,00 


rOO 


(1U5) 


-  ^  ^  ^2na   I       -I-    ^^.  „ia  [       ^L(^)0(.^).  £  (0(^).  |^(^(.^).  i) 


^-00 
,00 


v^Q 


'-00 


+  ka  Im  <       d^e 

,00 


2J4a 


.Ue--[iM^(^-|) 


(   .<=r-^  "  "^y  /.r-^       ..T 


V^I?      '^^        V^<2       k^ 


/2     2 
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At  this  juncture,  we  have  to  introduce  explicit  forms  of  9^,  and  it  is  first 
of  all  appropriate  to  utilize  the  approximation  7*  corresponding  to  infinite 
'aperture'  domain.  Now 


v^  v^W 


is  real  and  therefore  does  not  contribute  to  the  cross-section.  The  same  conclu- 
sion applies  to  the  integral 


r 

-00 


d^ 


Lv^ 


i>{xmO'  ^  (0fe)-  iVizfco*  0 


j       d^    1_ 


1- 


L   v^    /T 


because  its  branch  points  t-nd  poLes  at  ^  =   +  k  are  such  that  on  deforrdng  the 
path  of  integration  in  either  half -plane,    a  real  value  results.     In  its  simplified 
version,   then,    (lUS)  appears  as 


<r  «  na"^ 


/^     -     sin(2ka-  -A -k        -,        r  f    a         2ikaa    ,  > — ^_ 


yT? 


(Ili6) 


.00 


4  -    vC?^ 


with  an  obvious  change  of  variable  from  ^  to  x» 
the  latter  integrals,  for  ka  »  1,  yields 


An  asymptotic  evaluation  of 


(lh7) 


<r  it  na  |l 


L-2[v?-lj 


sin(2ka-  J) 


ka  »  1, 


-no- 


where the  numerical  coefficient,  2(  \fl  -l)  -  0,83,  of  the  second  term  is  to  be 
compared  with  the  values  1  and  2  obtained  by  the  refined  Kirchhoff  and  screen 
approaches  respectively.  The  'finite'  transform  |^(^),  (lUl),  obtained  by  inte- 
grating JTCx),  namely 


(^     .,                      ,    inA    ,     /^-^^^  it     -i^a+inA  /°°  it 
(1U8)      e-^^  ?J(x)dx  -  %  -  2_i L_  (      L_  dt-  t L.   e_  dt, 

k 

can  be  employed  in  the  cross-section  calculation,  but  this  does  not  reduce  the 
numerical  discrepancy.  It  would  seem  that  since  the  aperture  distribution  is 
accurate  near  the  rim  only,  the  dimension  of  the  aperture  is  effectively  accounted 
for  by  the  phase  factors  e~   ^  of  (1U5)« 

When  one  seeks  the  result  of  including  in  the  aperture  distribution  a  part 
which  arises  from  plane  wave  excitation  (e    ),  it  is  the  variational  formulation 
which  impedes  progress;  this  stems  from  the  fact  that  the  class  of  admissible  trial 
functions  }^  or  ^  is  required  to  take  a  given  value  at  the  rim.  Some  modification 
of  the  formulation  may  prove  advantageoxis,  but  will  not  be  considered  here, 

A  few  concluding  remarks  about  the  integral  equation  (llU)  are  appropriate. 
This  equation  may  be  handled  trivially  in  the  approximation  whereby  the  radical 

/2     2 

yk  -^     is  replaced  by  kj  the  Green's  fxinction  is  then  a  delta  fimction,   and  the 

solution  f  which  vanishes  at  p  ■  a  and  remains  finite  at  p  "  0  is  (cf.  (31)) 

J  (kp) 
(1U9)  I(p)  ■  1  -  f^  , 

from  (109),  we  deduce  the  corresponding  cross-section: 

(150)  Cr-  na2  -  -  ^^^  u    „a2(i  -  ^  tan(ka-  J)),      ka  »  1. 
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Note  that  the  wsve  function  (1U9)  tends  to  zero  in  a  linear  fashion  near  the  rim, 
while  at  the  aperture  center  the  imperturbed  part  is  generally  masked  by  large 
oscillations  of  the  remainder,  when  ka  »  1,  Furthemore,  if  for  the  inverse 
radical  we  write 


CtN'O  *^' 


(iDi)  can  be  converted  to  an  integral  equation  of  the  second  kind,  capable  of 
iteration.  It  appears,  however,  that  successive  iterations  have  some  effect  on 
the  first  correction  to  the  geometrical  cross-section,  and  hence  the  scheme  is 
not  really  practical. 


Acknowledgments 
I  am  greatly  indebted  to  Professors  Keller,  Marcuvitz  and  Schwinger  for 
valuable  discussions.  A  part  of  this  work  was  completed  during  tenure  of  a 
fellowship  of  the  John  Simon  Guggenheim  Memorial  Foundation. 


-  52  - 

Appendix  1 

In  Section  3j  the  screen  integral  equation  (25),  namely 


1-^  -^  I 
ik|p-p'| 


ik|p-p'|  _, 

^  .^^  .   dSt,       p  in  S 


_r,l 


(AI)  K(p')  dS«  . 

^S        |P-P'I        -^A   IP"P 

has  been  considered,  with  the  help  of  an  integral  representation  (28)  for  its 
kernel  (Green's  function).  Here  we  wish  to  indicate  another  approach  to  the 
integral  equation  and  cross-section  calculation.  Again,  the  circular  domain 
A  is  ass  lined,  and  if  p,  0  are  plane  polar  coordinates  with  origin  at  the  center, 
then  K  =  K(p),  while 


(A2)  \t-P\   -  /p^+p'^-2pp'  cos(0-^')  -  /(p-p')^+Upp'  sin^  2^  . 

The  straight- edge  approximation  (U2)  to  the  integral  equation  (Al)  follows  from 
the  restriction  |(i^'-/|  «  1,  which  implies  a  local  character  of  the  screen  dis- 
tribution. If  we  write  4  -  jf»  -  ^  and  set 


.  2  0'  «  ^   1  ^2 
sin   2^^  »  jj.  4 

the  angle  integrals  in  (Al)  become 


f  ik/(p-p')2+pp'c2 

j  dg, 

kl<A     /(p-p')^+  pp'   £,^ 
and  for  k-*  00,   this  can  be  estimated  by  the  corresponding  infinite  integral 


00 

(A3)  f 


exp[ik/(p-p.)2  +  pp.423 

— 'dK  »  — S^     Hi-'^(klp-p'l)   . 


/(p-p')2   *pp.  ^2  >/^P~ 
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Hence,    (Al)   takes  the  form 


00  a 

(Ah)  J     /F  K(p')H^^^k|p-p»|)dp'   -   f     y^H^^^klp-p'Ddp', 


p  >  a 


or,    using  the  integral  representation 


(A5)  H^^^k|p-p.|)  =i  j 


00 


-co  ykV 


we  have 


(A6) 


,00  ,00  ,00 

d^ 


'-00    yk"-^    ''a 


^f    V^K(p.)e^(p-P'^dp.  -f       -iL.e^pf)^e-^P'dp',        P  >  a. 


The  inhoraogeneous  term  in  (A6)  requires  further  approximation  to  secure  agreement 
with  Eq,  (U2),  and  we  note  first  that  it  can  be  written  equally  well  as 


.00 


/  ^\-^-\\     e^^M  ^    e-^P'dp., 


since  the  domains  of  p  and  p'  are  distinct.  Now,  the  significant  values  of  ^  are 
of  order  k,  and  asymptotic  development  for  the  p'  integral  yields 


j  VP-  e-^P'  dp.  -i^  I   ^— -dp 


=-i^P'  H.,  .  ^  d_  {     e-^P'  ... 


.  ■  d  r  V^  e-^"A   1  e-^%  1   ^e-^P*  .  .1 

^^L^? — ^— ^^i  ^;772H 


.i^e-.i^.O(A). 
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Thus,   retaining  only  the  leading  term,    (A6)  becomes 


(A7) 


^00  ,00 


.00       r 


•00 


v<?7 


/a  ^-00 


^VK? 


1 
1^ 


sin^(p-a),        p  >  a. 


when  the  odd  part  of  the  integrand  on  the  right-hand  side  of  (Ao)  is  dropped.  The 
change  of  variable  p  ■  a  •♦•  x,  and  use  of  (39),  (UO)  bring  (A?)  to  the  form  (U2). 

Next,  we  take  up  the  cross-section  calculation,  with  the  aim  of  utilizing 
the  straight- edge  distribution  (A7),  in  the  variational  form 


(A8)  Cr-G-j^+I^Im 


ik|p-p'| 


P-P'l 


K(P)  "        K(p')dSdS'-2    ^ 
u4       IP-P*I  ^s    IP-P'I 


dS« 


To  evaluate  the  first  of  two  angle  integrals  (ji^,0')  of  the  Green's  function,  we 
now  write 


(A9) 


^        expfik  /(p-p')^+Upp'sin^  2^ 


/(p-p')2  +  Upp-  sin^  £^ 


d^   +   J   d0 


exp 


ik/(p-p»)^  +  Upp'  sin^  2^1 
■      I  ■ 


v/(p-p»)^  +  Upp'  sin^  2^ 


and  adopt  the  previous  estimates  for  the  range  !0-0'|  <  A,  namely,  .     H^  (k|p-p'|X 

Vpp»" 

Since  in  the  remaining  interval  |^0*|  >  A  the  trigoncmetric  function  is  bounded 
away  from  zero,  and  since  k  is  large  there^the  radical  in  this  interval  is  approx- 
imated by  2  /pp'  I  sin  ^Ji    |  e 


-S$' 


Also 


\i 


f 


expfaik  v/ppT  I  sin  ||J 
2>/^  |sin|| 


exp 


2    I   d4 
o 

n 
2  j  d4 
o 


2ik  v^  sin  | 


exp 


ik 


2  Vpp^  sin  I 


exp 


[2ikv/^ 


sin 


d4 
exp 


VpF  k\ 


fpp"  i 


Ik/^  i 


2    ypp  •   sin 


ypp^  4 


*  0(f)' 


The  function 


■/: 


P(a)  -  I     cl4 
'o 


exp(2iasin  *)        io^ 

c  6 


2  sin 


can  be  shown  to  have  the  asymptotic  form 


P(a)^  l/f  exp[l(2a.«)]  .  0(i), 


and  consequently 


(AlO) 


^  exp[2ik  y^  I  sin  |l]  ^     expk(2k  y/^  -  J) 

1^-^^^  2^^^   |sin||  ^  (pp'K/'^ 
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Thus 


(All) 


.    I^a.  00 

K(p)  e"-^^^'^',  K(p')dSdS'  ^  ni(2n)   [  K(p)H^^^k|p-p'l  K(p' )  VpP'  dpdp' 

*  2n/|    e-"A|'^pp.)lAK(p)K(p.)e2i^^l^dpdp. 
/a 

00  /** 

^  2nia  [        -^^  I(^)I«)  *  2n  /f  expri(2ka-  J)l    (|    l(x)e^^dx)2 


,00 


2,^a  I       — S-  Ife)l(-i:)*  2«/¥  exp  i(2ka.  ?)     (I(-k))S 


where  we  have  approximated  exp  2ik  /ppM   by  exp  2ika(l+  )L   and  also  the 

geometrical  factor  1  +  —  by  unity.     Similsirly, 


(A12) 


K(p)dS   I :— dS'    ^ 


-2na 


4     lp-p'l 


r^^-2i"|/|e4(2.a.J)]l(-k). 


If  we  employ  the  transform  (of.    (5l)) 


i(^^(^)  -|(i-/rr|), 
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then 


cr^  <r^  .  ic/f  Jexp[i(2ka-  f)]  Jl<°)(-k)2  -  f  l'-°h-J 


■    ^K-|/f     =in(2ka.f) 


2i  na' 


X-i 


sin(2ka  -  r-X 
372—) 


Vn         (ka) 


ka  »  i 


which  is  the  well-known  result. 


Appendix  2 

We  consider  the  integral  equation 

00 

(A13)     [f(x»)H^-'-^(k|x-x'|)dx'  -  e"^^. 


X  >  0 


which  may  be  identified  as  the  equation  determining  the  discontinuity  in  normal 
derivative  of  the  wave  function,  f (x),  across  a  semi-infinite  soft  screen  occupying 
the  positive  x-axis,  excited  by  a  grazing  plane  wave  directed  at  the  straight  edge 
from  the  same  side  of  the  origin.  Eq,  (A13)  is  of  the  V/iener-Hopf  type,  and  pre- 


liminary to  solving  it,  we  write  an  extended  version: 


-ikx-6x 


00 

(Al)4)         j   f(x')H^'^^k|x-x'|)dx«  -< 

-00 


X  >  0 


'  ♦  g(x). 


X  <  0 
J 
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where  it  is  assumed  that 


(A15) 
and 


f  (x)  -  0, 


X  <  0 


(A16) 


g(x)  -  0, 


-6x 


X  >  0; 


the  convergence  factor,  e   ,  5  >  0,  has  been  introduced  in  the  plane  wave  term, 

and  its  relevance  will  be  seen  shortly.  If  both  sides  of  this  equation  are 

-i^x 
multiplied  by  e    ,  and  integrated  over  all  x,  it  follows  with  the  help  of  the 

relation 


that 


(A17) 


00 

H^l^k|x-x.|).l     f 

-oo 


^i^(x-x') 


<^* 


f(^) 


-^^*e(^^< 


Im  C 


-k+  iS 


k 

— •- 


ReC 


Figure  A 
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Now,  ifljnk"  e>0,  the   transform  of  the   Hankel  function  is  regular  in  the   strip 
|Iin  ^1  <  e,   the  transforms  of  f(x)  and  e  are  regular  in  the   (lower)  half -plane 

Im  K  <  -e+6,  while  the  transform  of  g(x)  is  regular  in  the   (upper)  half-plane  Im^  >  -c, 
The  common  region  of  regularity  defined  thereby,  namely  -e+6  >  Im  ^  >  -e,   owes  its 
finite  %ri.dth  to  the  parameter  6  (see  Fig.  U);  with  6-0,   the  upper  and  lower  half- 
planes  would  have  only  a  common  boundary  with  a  singular  point  on  it,  which  would 
make  for  complications.     Rewriting  the  transform  equation  as 

(A18)        ^*IVn    ^..\^^§^^\V^,iK), 

we  can  verify  that  the  left- and  right-hand  sides  are  regular  in  the  respective 
half-planes,  so  that  an  entire  function  is  thereby  defined  throughout  the  ^-plane. 
The  entire  function  actually  vanishes,  and  therefore 

(A19)  f  (^)  -  -  I   716  ^|L  . 

It  is  clear  that  both  the  transform  f(x)   and  its  inverse,  f (x),  tend  to  zero  with 
5.  However,  in  order  to  avoid  an  inconsistency  with  the  original  integral  equation, 
we  retain  6  /  0  and  first  calculate 


i 


00  ,y 

f(x')H^^\klx-x.|)dx.  -  I  I  eJ^lM 

-co  ^-00    ylc-c 


°°    i^ 


*-00 


where  the  latter  integration  path  lies  within  the  strip  of  regularity.  For  x  >  0, 
deformation  in  the  upper  half-plane  captures  the  pole  at  ^  ■  -k+i5,  and  the  integral 
is  simply  e      ,  which  reverts  to  an  undamped  plane  wave  if  6  ->  0  and  Im  k  ■  0. 
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